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Abstract 

In this paper, we exhibit a strong relation between the sand automata 
configuration space and the cellular automata configuration space. This 
relation induces a compact topology for sand automata, and a new context 
in which sand automata are homeomorphic to cellular automata acting 
on a specific subshift. We show that the existing topological results for 
sand automata, including the Hedlund-like representation theorem, still 
hold. In this context, we give a characterization of the cellular automata 
which are sand automata, and study some dynamical behaviors such as 
equicontinuity. Furthermore, we deal with the nilpotency. We show that 
the classical definition is not meaningful for sand automata. Then, we 
introduce a suitable new notion of nilpotency for sand automata. Finally, 
we prove that this simple dynamical behavior is undecidable. 



Keywords: sand automata, cellular automata, dynamical systems, subshifts, 
nilpotency, undecidability 

1 Introduction 

Self-organized criticality (SOC) is a common phenomenon observed in a huge 
variety of processes in physics, biology and computer science. A SOC system 
evolves to a "critical state" after some finite transient. Any perturbation, no 
matter how small, of the critical state generates a deep reorganization of the 
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whole system. Then, after some other finite transient, the system reaches a 
new critical state and so on. Examples of SOC systems are: sandpiles, snow 
avalanches, star clusters in the outer space, earthquakes, forest fires, load bal- 
ance in operating systems [21 13 HI 13 ED] . Among them, sandpiles models are a 
paradigmatic formal model for SOC systems [TT|[T2]. 

In [6] , the authors introduced sand automata as a generalization of sandpiles 
models and transposed them in the setting of discrete dynamical systems. A 
key-point of [6] was to introduce a (locally compact) metric topology to study 
the dynamical behavior of sand automata. A first and important result was a 
fundamental representation theorem similar to the well-known Hedlund's theo- 
rem for cellular automata [T3|. [6] . In [8], the authors investigate sand automata 
by dealing with some basic set properties and decidability issues. 

In this paper we continue the study of sand automata. First of all, we in- 
troduce a different metric on configurations (i.e. spatial distributions of sand 
grains) . This metric is defined by means of the relation between sand automata 
and cellular automata [8]. With the induced topology, the configuration set 
turns out to be a compact (and not only locally compact), perfect and to- 
tally disconnected space. The "strict" compactness gives a better topological 
background to study the behavior of sand automata (and in general of dis- 
crete dynamical systems). In fact, compactness provides a lots of very useful 
results which help in the investigation of several dynamical properties [H [16] . 
We show that all the topological results from [6] still hold, in particular the 
Hedlund-like representation theorem remains valid with the compact topology. 
Moreover, with this topology, any sand automaton is homeomorphic to a cel- 
lular automaton defined on a subset of its usual domain. We prove that it is 
possible to decide whether a given cellular automaton is in fact a sand automa- 
ton. Besides, this relation helps to prove some properties about the dynamical 
behavior of sand automata, such as the equivalence between equicontinuity and 
ultimate periodicity. 

Then, we study nilpotency of sand automata. The classical definition of 
nilpotency for cellular automata [TUHH] is not meaningful, since it prevents any 
sand automaton from being nilpotent. Therefore, we introduce a new definition 
which captures the intuitive idea that a nilpotent automaton destroys all the 
configurations: a sand automaton is nilpotent if all configurations get closer and 
closer to a uniform configuration, not necessarily reaching it. Finally, we prove 
that this behavior is undecidable. 

The paper is structured as follows. First, in Section [21 we recall basic defi- 
nitions and results about cellular automata and sand automata. Then, in Sec- 
tion [3l we define a compact topology and we prove some topological results, in 
particular the representation theorem. Finally, in Section^ nilpotency for sand 
automata is defined and proved undecidable. 
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2 Definitions 



For all a, b G Z with a < b, let [a, 6] = {a, a + 1, . . . , b} and [a, 6] = [a,b] U 
{+oo, — oo}. For a G Z, let [a, +oo) = {a, a + 1, . . .} \ {+oo}. Let N+ be the set 
of positive integers. For a vector i G Z d , denote by \i\ the infinite norm of i. 

Let A a (possibly infinite) alphabet and (f 6 N*. Denote by M d the set of 
all the d-dimensional matrices with values in A. We assume that the entries 
of any matrix U G M. d are all the integer vectors of a suitable d-dimensional 
hyper- rectangle [1, hi] X • • • X [1, hd] C N+. For any h= (hi,. . . ,h d ) eN|, let 
Mf, C M d be the set of all the matrices with entries in [1, hi] X • • • X [1, hd]. In 
the sequel, the vector h will be called the order of the matrices belonging to Mf. 
For a given element x £ A z , the /miie portion of x of reference position i G Z d 
and order h G is the matrix M£(x) G defined as Vfc G [1, /ii] x • • • x [1, h d ], 
MJ , l (x)k = Xi+k-i- For any r G N, let r d (or simply r if the dimension is not 
ambiguous) be the vector (r, . . . , r). 

2.1 Cellular automata and subshifts 

Let A be a finite alphabet. A CA configuration of dimension d is a function from 
Z d to A. The set A z of all the CA configurations is called the CA configuration 
space. This space is usually equipped with the Tychonoff metric At defined by 

Mx,y G A z \ d T (x,y)=2~ k where k = min{\j\ : j e Z d ,Xj ^ yj} . 

The topology induced by dr coincides with the product topology induced by 
the discrete topology on A. With this topology, the CA configuration space is a 
Cantor space: it is compact, perfect (i.e., it has no isolated points) and totally 
disconnected. 

For any k G Z d the shift map a k : A zd -> A zd is defined by Vx G A 1 ^ , Mi G 
Z d , a k (x)i = Xi+k- A function F : A z — > A z is said to be shift- commuting if 
Vfc GZ d , Foa k =a k oF. 

A d-dimensional subshift, S is a closed subset of the CA configuration space 
A zd which is shift-invariant, i.e. for any k G Z d , cr k (S) C S. Let T C and let 

be the set of configurations x € A 1, such that all possible finite portions of x 
do not belong to J 7 , i.e. for any i, h G Z d , M^(x) ^ JF. The set S 1 ^- is a subshift, 
and ^ is called its set of forbidden patterns. Note that for any subshift S, it is 
possible to find a set of forbidden patterns T such that S — Sp. A subshift S is 
said to be a subshift of finite type (SFT) if S = for some finite set T. The lan- 
guage of a subshift S is C(S) = {U G M d : 3i G Z d , h G N£, x g S*, M,^(x) = [/} 
(for more on subshifts, see [17] for instance). 

A cellular automaton is a quadruple (A, d, r, g) , where ^4 is the alphabet also 
called the state set, d is the dimension, r G N is the radius and g : M-2 r+ i —* A 
is the local rule of the automaton. The local rule g induces a global rule G : 
A 1 — ► A z defined as follows, 

Vx G A zd , Mi G 1 d , G(x)i = g{M£ r +1 (x)) . 
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Note that CA are exactly the class of all shift-commuting functions which are 
(uniformly) continuous with respect to the Tychonoff metric (Hedlund's theorem 
from [13]). For the sake of simplicity, we will make no distinction between a CA 
and its global rule G. 

The local rule g can be extended naturally to all finite matrices in the fol- 
lowing way. With a little abuse of notation, for any h £ [2r + 1, +oo) d and any 
U G Mf, define g(U) as the matrix obtained by the simultaneous application 
of g to all the A^r+i submatrices of U. Formally, g(U) = MJ i _ 2r (G(x)), where 
x is any configuration such that M^(x) = U. 

For a given CA, a state s G A is quiescent (resp., spreading) if for all matrices 
U G Mi r+1 such that Vfc G [1, 2r+l] d , (resp., 3k G [1, 2r+l] d ) U k = s, it holds 
that g(U) = s. Remark that a spreading state is also quiescent. A CA is said 
to be spreading if it has a spreading state. In the sequel, we will assume that 
for every spreading CA the spreading state is G A. 

2.2 SA Configurations 

A SA configuration (or simply configuration) is a set of sand grains organized in 
piles and distributed all over the d-dimensional lattice Z d . A pile is represented 
either by an integer from Z (number of grains), or by the value +00 (source 
of grains), or by the value — 00 (sink of grains), i.e. it is an element of Z = 
ZU{— 00, +00}. One pile is positioned in each point of the lattice Z d . Formally, 
a configuration a; is a function from Z^ to Z which associates any vector i = 
(ii, . . . , id) G Z d with the number Xi G Z of grains in the pile of position i. When 
the dimension d is known without ambiguity we note the null vector of 1 d . 
Denote by C = Z z the set of all configurations. A configuration x G C is said to 
be constant if there is an integer c G Z such that for any vector i G Z d , X{ = c. 
In that case we write x = c. A configuration x G C is said to be bounded if there 
exist two integers 7711,7712 G Z such that for all vectors i G Z d , mi < Xi < m%. 
Denote by B the set of all bounded configurations. 

A measuring device /3™ of precision r G N and reference height m G Z is a 
function from Z to [— r, r] defined as follows 



A measuring device is used to evaluate the relative height of two piles, with 
a bounded precision. This is the technical basis of the definition of cylinders, 
distances and ranges which are used all along this article. 

In [6], the authors equipped C with a metric in such a way that two con- 
figurations are at small distance if they have the same number of grains in a 
finite neighborhood of the pile indexed by the null vector. The neighborhood is 
individuated by putting the measuring device at the top of the pile, if this latter 
contains a finite number of grains. Otherwise the measuring device is put at 




if n > m + r , 
if 71 < 771 — r , 
otherwise. 
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height 0. In order to formalize this distance, the authors introduced the notion 
of cylinder, that we rename top cylinder. For any configuration x £ C, for any 
r £ N, and for any i £ Z d , the top cylinder of x centered in i and of radius r 
is the d-dimensional matrix C n r (x) £ M^ r+1 defined on the infinite alphabet 

A = Z by 

!Xi if k = r + 1 , 

0? (:r i+fc _ r _i) if fc ^ r + 1 and x, ^ ±oo , 
/3°(x i+ k_ r _i) otherwise. 

In dimension 1 and for a configuration x £ C, we have 

C r (x) = (/3 r 1 {xi— r ) , . . . , /? r 1 (x^_i ) , x 2 ; , f3 r 1 (x^-i-i ) , . . . , f3 r 1 (x^ r )) 
if ^ ±oo, while 

C' r (x) = (/3 r (Xj_ r ), ■ ■ ■ j P r {%i— l)j Xj, /3 r (Xj+l), • ■ ■ , /3,-(Xi+ r )) 

if x, = ±oo. 

By means of top cylinders, the distance d' : C x C — > M + has been introduced 
as follows: 

Vx,y£C, d'{x,y) = 2- k where fc = min {r e N : C'° r (x) ^ C'°(y)} . 

Proposition 2.1 (pH [8]) With the topology induced by d', the configuration 
space is locally compact, perfect and totally disconnected. 

2.3 Sand automata 

For any integer r £ N, for any configuration x £ C and any index i £ Z d 
with Xi ^ ±oo, the range of center i and radius r is the d-dimensional matrix 
r+1 on the finite alphabet A — [— r , r] U ± such that 

V* e[ l,2r + l] d , (*(*))* = { ^wist 1 ' 

The range is used to define a sand automaton. It is a kind of top cylinder, where 
the observer is always located on the top of the pile x, (called the reference). 
It represents what the automaton is able to see at position i. Sometimes the 
central ± symbol may be omitted for simplicity sake. The set of all possible 
ranges of radius r, in dimension d, is denoted by IZf. 

A sand automaton (SA) is a deterministic finite automaton working on con- 
figurations. Each pile is updated synchronously, according to a local rule which 
computes the variation of the pile by means of the range. Formally, a SA is a 
triple (d, r, f), where d is the dimension, r is the radius and / : IZf — > [—r, r] is 
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the local rule of the automaton. By means of the local rule, one can define the 
global rule F : C — » C as follows 

V.t e C, V^ e Z d , Fix), - > Xl 



x i + f{R l r { x )) otherwise. 

Remark that the radius r of the automaton has three different meanings: it rep- 
resents at the same time the number of measuring devices in every dimension of 
the range (number of piles in the neighborhood) , the precision of the measuring 
devices in the range, and the highest return value of the local rule (variation of 
a pile) . It guarantees that there are only a finite number of ranges and return 
values, so that the local rule has finite description. 

The following example illustrates a sand automaton whose behavior will be 
studied in Section [4j For more examples, we refer to [8]. 

Example 1 [the automaton N ] This automaton destroys a configuration by 
collapsing all piles towards the lowest one. It decreases a pile when there is a 
lower pile in the neighborhood (see Figure [J). Let Af = (1,1, /jy) of global rule 
Ftf where 

-1 if a < or b < , 



Va, b e [-1,1], fss(a,b) = 



otherwise. 

□ 
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Figure 1: Illustration of the behavior of J\f. 



When no misunderstanding is possible, we identify a SA with its global rule F. 
For any k & Z d , we extend the definition of the shift map to C, a k : C — > C 
is defined by Vx £ C, Vi G Z d , <J k (x)i = Xi + k- The raising map p : C — > C is 
defined by Vie e C, Vi G Z d , p(x)i = xi + 1. A function F : C — > C is said to be 
vertical-commuting \{ F o p = po F. A function F : C — > C is infinity-preserving 
if for any configuration idC and any vector i £ Z d , F(x)i = +oo if and only 
if = +oo and F{x)i = -co if and only if xi = -co. 

Remark that the raising map p is the sand automaton of radius 1 whose local 
rule always returns 1. On the opposite, the horizontal shifts cr, are not sand 
automata: they destroy infinite piles by moving them, which is not permitted 
by the definition of the global rule. 

Theorem 2.1 ([6, 8j) The class of SA is exactly the class of shift and vertical- 
commuting, infinity-preserving functions F : C — > C which are continuous w.r.t. 
the metric d'. 
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3 Topology and dynamics 



In this section we introduce a compact topology on the SA configuration space 
by means of a relation between SA and CA. With this topology, a Hedlund- 
like theorem still holds and each SA turns out to be homeomorphic to a CA 
acting on a specific subshift. We also characterize CA whose action on this 
subshift represents a SA. Finally, we prove that equicontinuity is equivalent to 
ultimate periodicity, and that expansivity is a very strong notion: there exist 
no positively expansive SA. 

3.1 A compact topology for SA configurations 

From [8], we know that any SA of dimension d can be simulated by a suitable CA 
of dimension d+1 (and also any CA can be simulated by a SA). In particular, 
a d-dimensional SA configuration can be seen as a (d + l)-dimensional CA 
configuration on the alphabet A = {0, 1}. More precisely, consider the function 

£ : C — > {0, 1} defined as follows 

A SA configuration x 6 C is coded by the CA configuration ((x) £ {0, 1} 
Remark that £ is an injective function. 

Consider the (d + l)-dimensional matrix K g M.f + 1, . . . , 1,2) such that 
-Ki,...,i,2 = 1 and = 0. With a little abuse of notation, denote Sk — 

S{k} the subshift of configurations that do not contain the pattern K. 

Proposition 3.1 The set ((C) is the subshift Sk ■ 

Proof. Each d-dimensional SA configuration x e C is coded by the (d + 1)- 
dimensional CA configuration ((x) such that for any i, h <E Af^(C(x)) ^ K, 

then £(C) C Sk- Conversely, we can define a preimage by ( for any y € Sk, by 
Vi G V 1 ,Xi = sup{fc : yu k) — !}• Hence ((C) = Sk- O 

Figure[2]illustrates the mapping £ and the matrix K = ( J for the dimension 

d = 1. The set of SA configurations C = Z z can be seen as the subshift 

z 2 

Sk = C,(C) of the CA configurations set {0, 1} . 

Definition 3.1 The distance d : C x C — > is defined as follows: 

Vx,yeC, d(x,y)=d T {((x),((y)) . 

In other words, the (well defined) distance d between two configurations x,y E C 
is nothing but the Tychonoff distance between the configurations ((x),((y) m 
the subshift Sk- The corresponding metric topology is the {0, 1} Z product 
topology induced on Sk- 
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(a) Valid configuration. 



(b) Invalid configuration. 



Figure 2: The configuration from Figure 2(a) is valid, while the configuration 
from Figure 2(b) contains the forbidden matrix K: there is a "hole". 



Remark 1 Note that this topology does not coincide with the topology obtained 
as countable product of the discrete topology on Z. Indeed, for any i £ Z d , 
the i th projection iti : C — > Z defined by iri(x) = Xi is not continuous in any 
configuration x with Xi = ±00. However, it is continuous in all configurations 
x such that Xi £ Z, since Vfc £ Z,Vx, y £ C, conditions 7r^(x) = k and d{x,y) < 

2 -max(| 4 |,fc) imply that n ( y } = fc 

By definition of this topology, if one considers ( as a map from C onto Sk , 
C turns out to be an isometric homeomorphism between the metric spaces C 
(endowed with d) and Sk (endowed with dr)- As an immediate consequence, 
the following results hold. 

Proposition 3.2 The set C is a compact and totally disconnected space where 
the open balls are clopen (i.e. closed and open) sets. 

Proposition 3.3 The space C is perfect. 

Proof. Choose an arbitrary configuration 16C. For any n £ N, let / £ Z d such 
that \l\ = n. We build a configuration y £ C, equal to a; except at site I, defined 
as follows 

w ■ rj,d \ m 1 I 1 if Xi = , 

VjeZ \{l}, yj = Xj and W = | Q otherwige _ 

By DefinitionEU d(y,x) = 2- n . □ 
Consider now the following notion. 

Definition 3.2 (ground cylinder) For any configuration x £ C, for any r £ 

N, and for any i £ 7L d , the ground cylinder of x centered on i and of radius r 
is the d- dimensional matrix C l r (x) £ A^2 r +i defined by 



Vfc 6 [l,2r + if, (C t r (x)) k = ffixi+k-r-!) 



For example in dimension 1, 

C l r (x) = (/3^(xj_ r ), . . . , (3 r (xi), . . . , /3 r (xj_|_ r )) . 

Figure [3] illustrates top cylinders and ground cylinders in dimension 1. Re- 
mark that the content of the two kinds of cylinders is totally different. 




(a) Top cylinder centered on Xi = 4: (b) Ground cylinder, at height 0: 

C'i.(x) = (+1, -oo, -3, 4, -2, -2, +1). C*(x) = (+oo, -2, +1, +oo, +2, +2, +oo). 

Figure 3: Illustration of the two notions of cylinders on the same configuration, 
with radius 3, in dimension 1. 

From Definition 13.11 we obtain the following expression of distance d by 
means of ground cylinders. 

Remark 2 For any pair of configurations x,y S C, we have 

d(x, y) = 2~ k where k = min {r e N : C°(x) ^ C®{y)} . 

As a consequence, two configurations x, y are compared by putting boxes (the 
ground cylinders) at height around the corresponding piles indexed by 0. 
The integer k is the size of the smallest cylinders in which a difference appears 
between x and y. This way of calculating the distance d is similar to the one 
used for the distance d', with the difference that the measuring devices and 
the cylinders are now located at height 0. This is slightly less intuitive than 
the distance d', since it does not correspond to the definition of the local rule. 
However, this fact is not an issue all the more since the configuration space 
is compact and the representation theorem still holds with the new topology 
(Theorem [32]). 

3.2 SA as CA on a subshift 

Let (X, mi) and (Y, 777.2) be two metric spaces. Two functions Hi : X — > X, 
H2 : Y — > Y are (topologically) conjugated if there exists a homeomorphism 
77 : X — > Y such that H2 ° rj = 77 o Hi . 
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We are going to show that any SA is conjugated to some restriction of a CA. 
Let F a d-dimensional SA of radius r and local rule /. Let us define the (d+ 1)- 
dimensional CA G on the alphabet {0, 1}, with radius 2r and local rule g defined 
as follows (see [8] for more details). Let M e -W^+i be a matrix on the finite 
alphabet {0, 1} which does not contain the pattern K. If there is a j 6 [r + 1, 3r] 

such that M(2r+i,...,2r+i,j) = 1 an d M(2r+i 2r+i,j+i) = 0, then let R € IZf 

be the range taken from M of radius r centered on (2r + 1, . . . ,2r + l,i). See 
figured] for an illustration of this construction in dimension d = 1. 




f 

Figure 4: Construction of the local rule g of the CA from the local rule / of 
the SA, in dimension 1. A range R of radius r is associated to the matrix M of 
order 4r + 1. 



The new central value depends on the height j of the central column plus its 
variation. Therefore, define g(M) = 1 if j + f(R) > 0, g(M) = if j + f(R) < 0, 
or g(M) = M( 2r +i i ... j 2r+i) (central value unchanged) if there is no such j. 

The following diagram commutes: 



4 1< > (d 

Sk ► Sk 

G 

i.e. Go( = (of. As an immediate consequence, we have the following result. 

Proposition 3.4 Any d- dimensional SA F is topologically conjugated to a suit- 
able (d + 1) -dimensional CA G acting on Sk- 

Being a dynamical submodel, SA share properties with CA, some of which 
are proved below. However, many results which are true for CA are no longer 
true for SA; for instance, injectivity and bijectivity are not equivalent, as proved 
in [7]. Thus, SA deserve to be considered as a new model. 

Corollary 3.3 The global rule F : C — > C of a SA is uniformly continuous w.r.t 
distance d. 

Proof. Let G be the global rule of the CA which simulates the given SA. Since 
the diagram {J} commutes and £ is a homeomorphism, F = oGo(. Since 
G is a continuous map and, by Proposition 13.2) C is compact, then the thesis is 
obtained. □ 
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For every a G Z, let P a = 7r 1 ({ a l) be the clopen (and compact) set of all 
configurations x G C such that xq = a. 

Lemma 3.4 Let F : C — » C be a continuous and infinity-preserving map. There 
exists an integer I G N such that for any configuration x G Pq we have |.F(:r)o| < 

Proof. Since F is continuous and infinity-preserving, the set F(Pq) is compact 
and included in 7r ( ^" 1 (Z). From Remark [TJ ttq is continuous on the set 7r _1 (Z) 
and in particular it is continuous on the compact F(Pq). Hence ttq(F(Pq)) is a 
compact subset of Z containing no infinity, and therefore it is included in some 
interval [—1, I], where I G N. □ 

Theorem 3.5 A mapping F : C — > C is the global transition rule of a sand 
automaton if and only if all the following statements hold 

(i) F is (uniformly) continuous w.r.t the distance d; 
(ii) F is shift- commuting; 
(Hi) F is vertical-commuting; 
(iv) F is infinity-preserving. 

Proof. Let F be the global rule of a SA. By definition of SA, F is shift- 
commuting, vertical-commuting and infinity-preserving. From Corollary 13,31 F 
is also uniformly continuous. 

Conversely, let F be a continuous map which is shift-commuting, vertical- 
commuting, and infinity-preserving. By compactness of the space C, F is also 
uniformly continuous. Let I G N be the integer given by Lemma I3T41 Since F is 
uniformly continuous, there exists an integer r G N such that 

Vx, y G C C° r (x) = C° r (y) => C?(F(x)) - Cf(F(y)) . 

We now construct the local rule / : IZf — ► [— r, r) of the automaton. For any 
input range R G TZf, set f(R) = F(x)o, where x is an arbitrary configuration of 
P Q such that Vfc G [1, 2r + 1], k ^ r+1, (3°(x k - r -i) = R k . Note that the value of 
f(R) does not depend on the particular choice of the configuration x G Pq such 
that Vfc 7^ r+ 1, /3°(xk- r -i) = Rk- Indeed, Lemma [3~4l and uniform continuity 
together ensure that for any other configuration y G Pq such that Vfc ^ r + 1, 
P° r (y k - r -i) = Rk, we have F(y) - F(x) , since $(F(x) ) = P?(F(y) ) and 
|-F 1 (y)o| < l- Thus the rule / is well defined. 

We now show that F is the global mapping of the sand automaton of radius 
r and local rule /. Thanks to (JTvJ) , it is sufficient to prove that for any x G C 
and for any i G Z d with \x%\ ^ oo, we have F(x)i = Xi + f (R l r (x)) . By Jn| and 
Jm|, for any i G Z d such that \xi\ ^ oo, it holds that 

F(x) i =[p^oa- i (F(a i op-^(x)))] i 

= x l+ [a- l (F(o-*op-*>(x)))} % 

= x i +[F(a i op-^(x))] Q . 
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Since a 1 o p Xi (x) G Pq, we have by definition of / 

F(x) i = x i + f(R°(a i op- x *(x))) . 

Moreover, by definition of the range, for all k G [1, 2r + l] d , 

R° r ^op-^{x)) k = (5f°P^'^^op-^{x) k ) = o r {x i+k -Xi) = ft*(x i+1c ) , 

hence R° r {a l o p~ Xi (x)) = R^x), which leads to F{x) l =Xi + f (R l r {x)). □ 

We now deal with the following question: given a (d + l)-dimensional CA, 
does it represent a d-dimensional SA, in the sense of the conjugacy expressed 
by diagram Q]? In order to answer to this question we start to express the 
condition under which the action of a CA G can be restricted to a subshift Sjr, 
i.e., G(Sjr) C Sjr (if this fact holds, the subshift Sjf is said to be G-invariant). 

Lemma 3.6 Let G and Sp be a CA and a subshift of finite type, respectively. 
The condition G(Sp) C is satisfied iff for any U G C(S^) and any H G T 
of the same order than g{U), it holds that g(U) ^ H . 

Proof. Suppose that G(Sjr) C Syr. Choose arbitrarily H G T and U G C(Sj^), 
with g{U) and H of the same order. Let x G Sp containing the matrix U. Since 
G{x) G Sp, then g(U) G C(Sj^), and so g(U) ^ H. Conversely, if x G Sj? and 
G(x) i Sjr, then there exist U G £(S» and H G T with g(U) = H. □ 

The following proposition gives a sufficient and necessary condition under which 
the action of a CA G on configurations of the G-invariant subshift Sk = C 
preserves any column whose cells have the same value. 

Lemma 3.7 Let G be a (d+1) -dimensional CA with state set {0, 1} and Sk be 
the subshift representing SA configurations. The following two statements are 
equivalent: 

(i) for any x G Sk with X(o,....o.;) = 1 (resp., a5(o,...,o,») — ®) f or a ^ * G Z, it 
holds that G(cc)( ,....o.i) = 1 (resp., G(x)(o,...,o,i) = 0) for all i £ Z. 

(ii) for any matrix U G M^t+i ^ £{Sk) with f/( r +i,...,r+i,fe) = 1 (resp., 
U( r +i r +i k) = 0) o,nd any k G [1, 2r + 1], it holds that g{U) = 1 (resp., 
g(U) = 0). ' 

Proof. Suppose that (1) is true. Let U G M^ r+ i H C(Sk) be a matrix with 
^(r+i,...,r+i,fc) = 1 an d l e t x G Sk be a configuration such that £(o,...,o,i) = 1 
for all £ G Z and M^ r r ) _ 1 (a;) = {/. Since G(x)(p,...,o,i) = 1 for all i G Z, and 
^2r+i( a: ) — U, then <?(£/) = 1. Conversely, let x G Sk with £(o,...,o,i) = 1 for 
all i G Z. By shift-invariance, we obtain G(x)(o,...,o,i) = 1 for all i G Z. □ 

Lemmas 13.61 and 13771 immediately lead to the following conclusion. 

Proposition 3.5 It is decidable to check whether a given (d + 1) -dimensional 
CA corresponds to a d- dimensional SA. 
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3.3 Some dynamical behaviors 

SA are very interesting dynamical systems, which in some sense "lie" between 
d-dimensional and d + 1-dimensional CA. Indeed, we have seen in the previ- 
ous section that the latter can simulate d-dimensional SA, which can, in turn, 
simulate d-dimensional CA. For the dimension d = 1, a classification of CA in 
terms of their dynamical behavior was given in [15]. Things are very different 
as soon as we get into dimension d = 2, as noted in [19J [18] . The question is 
now whether the complexity of the SA model is closer to that of the lower or 
the higher-dimensional CA. 

Let (X, m) be a metric space and let H : X — > X be a continuous appli- 
cation. An element x £ X is an equicontinuity point for H if for any e > 0, 
there exists S > such that for all y £ X, m(x,y) < S implies that Vn £ N, 
m(H n (x) , if" (y)) < e. The map H is equicontinuous if for any e > 0, there 
exists 5 > such that for all x, y £ X, m(x,y) < S implies that Vn £ N, 
m(H n (x), H n (y)) < e. If X is compact, H is equicontinuous iff all elements of 
X are equicontinuity points. An element x £ X is ultimately periodic for H if 
there exist two integers n > (the preperiod) and p > (the period) such that 
H n+P (x) = H n (x). H is ultimately periodic if there exist n > and p > such 
that H n+P = H n . H is sensitive (to the initial conditions) if there is a constant 
e > such that for all points x € X and all 8 > 0, there is a point y £ X and an 
integer n £ N such that m{x,y) < 8 but m{F n (x) , F n (y)) > e. H is positively 
expansive if there is a constant e > such that for all distinct points x, y £ A, 
there exists rt £ N such that m(H n (x) , H n (y)) > e. 

The topological conjugacy between a SA and some CA acting on the special 
subshift Sk helps to adapt some properties of CA. In particular, the following 
characterization of equicontinuous C A can be adapted from Theorem 4 of [15] . 

Proposition 3.6 If F is a SA, then the following statements are equivalent: 

1. F is equicontinuous. 

2. F is ultimately periodic. 

3. All configurations of C are ultimately periodic for F. 

Proof. d^H For any n > and p > 0, let D n , p = {x : F n+P (x) = F n (x)}. 
Remark that C = \J n pgN D„ jP is the union of these closed subsets. As C is com- 
plete of nonempty interior, by the Baire Theorem, there are integers n,p £ N for 
which the set D„ jP has nonempty interior. Hence the conjugate image Q(D n ^ p ) 
has nonempty interior too, and it can easily be seen that it is a subshift. It 
is known that the only subshift with nonempty interior is the full space; hence 
D — C 
H^GD obvious. 

[2]=^lJ Let F be ultimately periodic with F n+P = F n for some n > 0, p > 0. 
Since F,F 2 ,... ; F n +P _1 are uniformly continuous maps, for any e > there 
exists S > such that for all x,y £ C with d(x, y) < S, it holds that \fq £ N, 
q < n+p, d(F q (x) , F q (y)) < e. Since for any t £ N F* is equal to some F q with 
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q < n + p, the map F is equicontinuous. 

[T]=^2J For the sake of simplicity, we give the proof for a given one-dimensional 
equicontinuous SA F. Let G be the global rule of the two-dimensional CA 
whose action on Sk is conjugated to F. By Definition 13. 1\ and since the 
diagram Q] commutes, the map G : Sk — > Sk is equicontinuous w.r.t. dy. 
So, for e = 1, there exists I G N such that for all x, y G Sjf, if M^ +1 (x) = 
M 2i+i(y)> then for a11 * e N > G*(a;) = G t (y) . Consider now configurations 
C(c), where c G {— oo, +oo} z has either the form (. . . , — oo, — oo, +00, +00, . . .) 
or (. . . , +00, +00, —00, — 00, . . .). Since every ((c) are ultimately periodic (with 
preperiod n = and period p = 1) and G is equicontinuous, for any k G Z 2 
and any y G Sr- with M 2 ^, 1 (y) = M 2 ^ +1 (((c)), it holds that the sequence 
{G*(y)fc}t S N is ultimately periodic. For any U G C(Sk) HM^j, let x u be 
the configuration such that M^ +1 (x) = U, = if — / < i < I and 

j > I, and xuj\ = 1 otherwise. Except for the finite central region, x u is 
made by the repetition of a finite number of matrices appearing inside con- 
figurations C( c )- Hence, x u is an ultimately periodic configuration with some 
preperiod n\j and period pjj. Then, for any y G Sk with M^ +1 (y) = U, the 
sequence {G*(y)o}teN is ultimately periodic with preperiod njj and period pu ■ 
Set n = max{n!7 : U G C(S K )r)Ml l+1 } andp = lcm{pu : U G C(S K )r\M\ l+1 \ 
where 1cm is the least common multiple. Thus, for any configuration z G Sk, 
we have that G n (z)o = G n+p (z)o. By shift-invariance, we obtain Vfc G Z 2 , 
G n (z)k = G n+P (z)k. Concluding, G is ultimately periodic and then F is too. □ 

In [15] is presented a classification of CA into four classes: equicontinuous 
CA, non equicontinuous CA admitting an equicontinuity configuration, sensitive 
but not positively expansive CA, positively expansive CA. This classification is 
no more relevant in the context of SA since the class of positively expansive 
SA is empty. This result can be related to the absence of positively expansive 
two-dimensional CA (see [3jj]), though the proof is much different. 

Proposition 3.7 There are no positively expansive SA. 

Proof. Let F a SA and 5 = 2~ k > 0. Take two distinct configurations x, y G C 
such that Vi G [—k,k],Xi = y t = +00. By infinity-preservingness, we get 
Vn G N,Vi G [-k,k],F n (x)i = F n (y) t = +00, hence d(F n (x) , F n (y)) < S. □ 

An important open question in the dynamical behavior of SA is the existence 
of non-sensitive SA without any equicontinuity configuration. An example for 
two-dimensional CA is given in [18], but their method can hardly be adapted 
for SA. This could lead to a classification of SA into four classes: equicon- 
tinuous, admitting an equicontinuity configuration (but not equicontinuous), 
non-sensitive without equicontinuity configurations, sensitive. 

Another issue is the decidability of these classes. In [7], the undecidability 
of SA ultimate periodicity was proved on the particular subsets of finite and 
periodic configurations. It follows directly that equicontinuity on these subsets 
is undecidable. The question is still open for the whole configuration space C. 
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4 The nilpotency problem 



In this section we give a definition of nilpotency for SA. Then, we prove that 
nilpotency behavior is undecidable f Theorem I4.5p . 

4.1 Nilpotency of CA 

Here we recall the basic definitions and properties of nilpotent CA. Nilpotency 
is among the simplest dynamical behavior that an automaton may exhibit. In- 
tuitively, an automaton defined by a local rule and working on configurations 
(either C or A 1, ) is nilpotent if it destroys every piece of information in any 
initial configuration, reaching a common constant configuration after a while. 
For CA, this is formalized as follows. 

Definition 4.1 (CA nilpotency |10L I14j) A CA G is nilpotent if 

3c e A, 3N 6 N Vie A zd , Vn > N, G n (x) = c . 

Remark that in a similar way to the proof of Proposition 13.61 Definition 14.11 
can be restated as follows: a CA is nilpotent if and only if it is nilpotent for all 
initial configurations. 

Spreading CA have the following stronger characterization. 

Proposition 4.1 (|9j) A CA G, with spreading state 0, is nilpotent iff for every 
x £ A z , there exists n £ N and i £ Z d such that G n (x)i =0 (i.e. appears in 
the evolution of every configuration) . 

The previous result immediately leads to the following equivalence. 

Corollary 4.2 A CA of global rule G, with spreading state 0, is nilpotent if 
and only if for all configurations x £ A z , linin^oo dr(G n (a;), 0) = 0. 

Recall that the CA nilpotency is undecidable [14]. Remark that the proof 
of this result also works for the restricted class of spreading CA. 

Theorem 4.3 ([14]) For a given state s, it is undecidable to know whether a 
cellular automaton with spreading state s is nilpotent. 

4.2 Nilpotency of SA 

A direct adaptation of Definition [4J] to SA is vain. Indeed, assume F is a SA of 
radius r. For any k £ Z d , consider the configuration x k £ B defined by Xq = k 
and x\ = for any i £ Z d \{0}. Since the pile of height k may decrease at most 
by r during one step of evolution of the SA, and the other piles may increase at 
most by r, x k requires at least \k/2r~\ steps to reach a constant configuration. 
Thus, there exists no common integer n such that all configurations x k reach 
a constant configuration in time n. This is a major difference with CA, which 
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is essentially due to the unbounded set of states and to the infinity-preserving 
property. 

Thus, we propose to label as nilpotent the SA which make every pile ap- 
proach a constant value, but not necessarily reaching it ultimately. This nilpo- 
tency notion, inspired by Proposition 14.21 is formalized as follows for a SA F: 

3c G Z, Vie C, lim d(F"(x),c) = . 

n — >oo 

Remark that c shall not be taken in the full state set Z, because allowing in- 
finite values for c would not correspond to the intuitive idea that a nilpotent 
SA "destroys" a configuration (otherwise, the raising map would be nilpotent). 
Anyway, this definition is not satisfying because of the vertical commutativity: 
two configurations which differ by a vertical shift reach two different configura- 
tions, and then no nilpotent SA may exist. A possible way to work around this 
issue is to make the limit configuration depend on the initial one: 

Vx G C, 3c e Z, lim d(F n (x),c) = . 

n — >oc 

Again, since SA are infinity-preserving, an infinite pile cannot be destroyed 
(nor, for the same reason, can an infinite pile be built from a finite one). There- 
fore nilpotency has to involve the configurations of Z z , i.e. the ones without 
infinite piles. Moreover, every configuration x G Z z made of regular steps (i.e. 
in dimension 1, for all i G Z, X{ — Xj_i = x,+i — Xi) is invariant by the SA rule 
(possibly composing it with the vertical shift). So it cannot reach nor approach 
a constant configuration. Thus, the larger reasonable set on which nilpotency 
might be defined is the set of bounded configurations B. This leads to the 
following formal definition of nilpotency for SA. 

Definition 4.4 (SA nilpotency) A SA F is nilpotent if and only if 
Vx G B, 3c G Z, lim d(F n (x),c) = . 

n — >oo 

The following proposition shows that the class of nilpotent SA is nonempty. 

Proposition 4.2 The SA Af from Example^ is nilpotent. 

Proof. Let x G B, let i G Z such that for all j G Z, Xj > X{. Clearly, after 
Xi+i — Xi steps, Fp +1 ~ Xz (x)i+i = F^ +1 ~ x ' (x)i = x,. By immediate induction, 
we obtain that for all j G Z there exists rij G N such that Fjj{x)j = Xi, hence 
lim„^ 00 d(F^(x),x i ) = 0. □ 

Similar nilpotent SA can be constructed with any radius and in any dimen- 
sion. 
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4.3 Undecidability 

The main result of this section is that SA nilpotency is undecidable (Theo- 
rem [431), by reducing the nilpotency of spreading CA to it. This emphasizes 
the fact that the dynamical behavior of SA is very difficult to predict. We think 
that this result might be used as the reference undecidable problem for further 
questions on SA. 

Problem Nil 

instance: a SA A — (d, r, A); 
QUESTION: is A nilpotent? 

Theorem 4.5 The problem Nil is undecidable. 

Proof. This is proved by reducing Nil to the nilpotency of spreading cellular 
automata. Remark that it is sufficient to show the result in dimension 1. Let S 
be a spreading cellular automaton S = (A, l,s,g) of global rule G, with finite 
set of integer states 4cN containing the spreading state 0. We simulate S with 
the sand automaton A = (l,r = max(2s, max A), f) of global rule F using the 
following technique, also developed in [8]. Let £ : A z — ► B be a function which 
inserts markers every two cells in the CA configuration to obtain a bounded 
SA configuration. These markers allow the local rule of the SA to know the 
absolute state of each pile and behave as the local rule of the CA. To simplify 
the proof, the markers are put at height (see Figured]): 

Vy e A\ V, e Z, Z(v)i = I ° (marker) * [ iS ° dd ' 
y ' ' \ y i/2 otherwise. 

This can lead to an ambiguity when all the states in the neighborhood of size 
4s + 1 are at state 0, as shown in the picture. But as in this special case the 
state is quiescent for g, this is not a problem: the state is preserved, and 
markers are preserved. 




f 

Figure 5: Illustration of the function £ used in the simulation of the spreading 
CA S by A. The thick segments are the markers used to distinguish the states 
of the CA, put at height 0. There is an ambiguity for the two piles indicated 
by the arrows: with a radius 2, the neighborhoods are the same, although one 
of the piles is a marker and the other the state 0. 
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The local rule / is defined as follows, for all ranges R G IZl, 



if -R_2 S +i,-R-2s+3) • • • ,R-i, Ri, ■ ■ • i -R2S-1 £ A , 

f(R) = { g{R-2s + a, R-2.S+2 + a, ■ ■ ■ , R-2 + a,a,R 2 + a, . . . , R 2s + a) - a 

if R-2s+i = R-2s+3 = ■ ■ ■ = R2.S-1 = a < and -a G A . 

(2) 

The first case is for the markers (and state 0) which remain unchanged, the 
second case is the simulation of g in the even piles. As proved in [8], for any 
y G A z it holds that £{G(y)) = F(£(y)). The images by / of the remaining 
ranges will be defined later on, first a few new notions need to be introduced. 

A sequence of consecutive piles (xj, . . . , Xj) from a configuration x G B is 
said to be valid if it is part of an encoding of a CA configuration, i.e. Xi = 
Xi + 2 = ■ ■ ■ = Xj (these piles are markers) and for all k G N such that < k < 
(j — i)/2, Xi+2k+i — Xi £ A (this is a valid state). We extend this definition to 
configurations, when i = —00 and j = +00, i.e. x € p c o £(A Z ) for a given c 6 Z 
(x G B is valid if it is the raised image of a C A configuration) . A sequence (or 
a configuration) in invalid if it is not valid. 

First we show that starting from a valid configuration, the SA A is nilpotent 
if and only if S is nilpotent. This is due to the fact that we chose to put the 
markers at height 0, hence for any valid encoding of the CA x = p c o £(y), with 
y e A z and c G Z, 

lim d T (G"(y),0) = if and only if lim d(F n {x), c) = . 

n — >oc n — >oc 

It remains to prove that for any invalid configuration, A is also nilpotent. In 
order to have this behavior, we add to the local rule / the rules of the nilpotent 
automaton Af for every invalid neighborhood of width 4s + 1. For all ranges 
R G TZl not considered in Equation ((2]), 

_ J —1 if R-r < or R-r+i < or ••■ or Rr < , , . 

^ H >~{ otherwise. (6 > 

Let x G B be an invalid configuration. Let k G Z be any index such that 
VZ G Z, xi > Xk- Let i,j £ Z be respectively the lowest and greatest indices 
such that i < k < j and (xj, . . . , Xj) is valid (i may equal j). Remark that 
for all n G N, (F"(a:)j, . . . ,F n (x)j) remains valid. Indeed, the markers are by 
construction the lowest piles and Equations ^ and (O do not modify them. 
The piles coding for non-zero states can change their state by Equation 
or decrease it by 1 by Equation ([3]), which in both cases is a valid encoding. 
Moreover, the piles Xi-i and Xj+i will reach a valid value after a finite number 
of steps: as long as they are invalid, they decrease by 1 until they reach a value 
which codes for a valid state. Hence, by induction, for any indices a, b G Z, 
there exists N a ^ such that for all n > N a ^ the sequence (F n (x) a , . . . , F n {x)b) 
is valid. 

In particular, after N_2Nr-i,2Nr+i step, there is a valid sequence of length 
ANr + 3 centered on the origin (here, N is the number of steps needed by S to 
reach the configuration 0, given by Definition 14. ip . Hence, after N-2Nr,2Nr + N 
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steps, the local rule of the CA S applied on this valid sequence leads to 3 
consecutive zeros at positions —1, 0, 1. All these steps are illustrated on Figure[6l 



i j 




















1 1 1 










1 1 : 







time 



Figure 6: Destruction of the invalid parts. The lowest valid sequence (in gray) 
extends until it is large enough. Then after N other steps the 3 central piles 
(hatched) are destroyed because the rule of the CA is applied correctly. 



Similarly, we prove that for all n > N-2Nr-k,2Nr+k + N, the sequence 
(F n (x)-k, ■ ■ ■ , F n (x)k) is a constant sequence which does not evolve. There- 
fore, there exists c G Z such that linin^oo d(F n (x),c) — 0. We just proved that 
A is nilpotent, i.e. \im n -, 00 d(F n (x),c) = for all x £ B, if and only if S is 
nilpotent (because of the equivalence of definitions given by Corollary I4.2j) , so 
Nil is undecidable (Proposition 14.3ft . □ 



5 Conclusion 

In this article we have continued the study of sand automata, by introducing a 
compact topology on the SA. In this new context of study, the characterization 
of SA functions of [6j [8| still holds. Moreover, a topological conjugacy of any 
SA with a suitable CA acting on a particular subshift might facilitate future 
studies about dynamical and topological properties of SA, as for the proof of the 
equivalence between equicontinuity and ultimate periodicity (Proposition I3.6p . 

Then, we have given a definition of nilpotency. Although it differs from the 
standard one for CA, it captures the intuitive idea that a nilpotent automaton 
"destroys" configurations. Even though nilpotent SA may not completely de- 
stroy the initial configuration, they flatten them progressively. Finally, we have 
proved that SA nilpotency is undecidable (Theorem I4.5|l . This fact enhances 
the idea that the behavior of a SA is hard to predict. We also think that this 
result might be used as a fundamental undecidability result, which could be 
reduced to other SA properties. 
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Among these, deciding dynamical behaviors remains a major problem. More- 
over, the study of global properties such as injectivity and surjectivity and 
their corresponding dimension-dependent decidability problems could help un- 
derstand if d-dimensional S A look more like d-dimensional or d + 1-dimensional 
CA. Still in that idea is the open problem of the dichotomy between sensitive 
SA and those with equicontinuous configurations. A potential counter-example 
would give a more precise idea of the dynamical behaviors represented by SA. 
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